PUPT-1503 
hep-th/9410054 



Some Properties of Large N 
^ ■ Two Dimensional Yang-Mills Theory 



o 
O 



> 

o 
o 



David J. Gross''' 



G^ ■ gross@puh.epl.princeton.edu 



and 

Andrei Matytsin 
matytsinOpuhepl . princeton . edu 



. Department of Physics 

Joseph Henry Laboratories 
p • Princeton University 

^ ; Princeton, NJ 08544 

Abstract 

^ • Large two-dimensional QCD on a cylinder and on a vertex manifold (a sphere with 
three holes) is investigated. The relation between the saddle-point description and the 
collective field theory of QCD2 is established. Using this relation, it is shown that the 
Douglas-Kazakov phase transition on a cylinder is associated with the presence of a gap 
in the eigenvalue distributions for Wilson loops. An exact formula for the phase transition 
on disc with an arbitrary boundary holonomy is found. The role of instantons in inducing 
such transitions is discussed. The zero-area limit of the partition function on a vertex 
manifold is studied. It is found that this partition function vanishes unless the boundary 
conditions satisfy a certain selection rule which is an analogue of momentum conservation 
in field theory. 
September 1994 



^ This work was supported in part by the National Science Foundation under grant PHY90- 
21984. 



1. Introduction. 

Recently there has been much interest in the study of QCD2 in the large limit, 
largely motivated by an attempt to find a string representation of QCD in four dimensions 
0]. The study of large QCD2 is also useful in exploring more genral properties of large 
QCD. As such, a detailed analysis of the properties of the theory on the sphere and on 
the cylinder are interesting. Perhaps the most important of these properties is the third 
order phase transition for large QCD on a sphere, discovered recently by Douglas and 
Kazakov p|. Physically, this transition is caused by condensation of instantons. One can 
demonstrate that the effect of instantons is negligible when the area of the sphere is small, 
and becomes significant as the area reaches the critical value, Acr = tt^/A, where A = g'^N 
is the large A^ coupling constant of QCD 0. 

The existence of this phase tranition, which is similar to the phase transition that 
occurs in the one-plaquette model 0, is of great interest. Such effects signal a sharp tran- 
sition between weak and strong coupling, between non-confining physics and confinement, 
and between the a stringy regime and a non-string regime. If they were to occur in four 
dimensional QCD they might be an indication that the string picture is limited to infrared 
phenomena and cannot be analytically continued to the ultraviolet domain. One of our 
goals is to analyse the occurence of such phase transitions for other geometries. We will 
develop methods that allow us to explore this phase transition in detail. 

A phase transition similar to the Douglas-Kazakov transition occurs in QCD2 on a 
cylinder with fixed boundary conditions (that is, with fixed holonomy matrices Uc^ = 
Vexp A^{x) dx^ and = "Pexp A^{x) dx^ . Here Ci and C2 are the two circles 
forming the boundary of the cylinder). Obviously, such a phase transition should manifest 
itself in the master field of two-dimensional QCD []5|. 

A natural set of physical observables in gauge theory is formed by the Wilson loops 
Wn{C) = iTr with Uc = Vexp §^ A^{x) dx^" . Therefore, we could say that the master 
field is fully described by the set of quantities Wn{C) for all possible contours C. Indeed, 
Wn{C) satisfy a closed set of equations (the loop equations) with a well defined large A^ 
limit In addition, the expectation values of Wilson loops do factorize at large A^. 

The loop equations, however, are higly nonlocal and the description of the theory they 
provide is very convoluted. 

The loop variables are difficult to work with. For two-dimensional QCD a better set 
of variables is readily available. Indeed, it is known that the two-dimensional QCD is 
equivalent to a c = 1 matrix model with the spatial coordinate compactified on a circle 
Therefore, the eigenvalue density of the Wilson matrix Uc will satisfy the 
appropriate collective field equation (the Hopf equation) as a function of the area bounded 
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by the contour C. Using this fact we can obtain a formula for Wilson loops Wn in the 
large N QCD on a cylinder with arbitrary boundary conditions. 

As a result, the dynamics of Wilson loops on a cylinder can be desribed by a simple 
physical picture. Let cri{6) and (72(6'), 9 G [0,27r], be the eigenvalue densities of the 

boundary holonomy matrices Uc^ and Uc2- Consider a one- dimensional compressible fluid 

2 „ 

living on a circle, with the equation of state P = — i^a , P being the pressure and a the 
density. Imagine the process where the fluid moves from an initial conflguration, where its 
density proflle is (j[6) = a'i{6) to a flnal one, (j{6) = 0-2(0) , during the time interval equal 
to the total area of the cylinder, A. Then the eigenvalue density for the Wilson matrix 
Uc, where the contour C cuts the cylinder of total area A into two cylinders of areas Ai 
and A2, is just given by the density of our fluid at the time Ai. Consequently, 

Wn{C) = J a{t^Ai,e) e'""^ dO. (1.1) 


On the other hand, the partition function of the two-dimensional QCD, at least in 
the most interesting case of symmetric boundary conditions Uc^ = Uq^^ is dominated by 
a single representation of the gauge group. This representation, associated with its Young 
tableau, can be characterized by a set of numbers hi = k/N, i = 1, . . . , A?". The density 
of this set, pvih), plays the crucial role in the analysis of the Douglas-Kazakov phase 
transition. However, there is no known equation for determining pvih), except for some 
special cases. 



Ai A2 




A/2 A/2 

Fig. 1: The cylinder with boundary contours Ci and C2. The contour Co, used 
in the calculation of the dominant Young tableau, cuts the cylinder in half. 
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This problem shall be solved below, in section 2. We find that, when Uc^ = Uq^, the 
Young tableau density pvih) is determined by the surprisingly simple formula 

npY{-nao{e))=e (1.2) 

where ao{9) is the density of eigenvalues for the Wilson matrix Ucg- Here the contour Co 
cuts the cylinder into two equal parts, Ai = A2 = A/2 (see fig. 1). The density ao{9) can 
be found using the procedure outlined aboveB. This formula will also allow us to prove 
the criterion (formulated by Caselle, D'Adda, Magnea and Panzeri |jTl|]) determining for 
which boundary conditions the partition function of QCD2 can exhibit a Douglas-Kazakov 
phase transition, and for which the transition is absent. This criterion applies even in those 
situations when no single representation is dominant and the original method of Douglas 
and Kazakov leads to a complex saddle point. Quite remarkably, we discover that for QCD 
on a disc (that is, when a2{9) = S{9)) the transition point can be determined exactly and 
explicitly to be 

ai{9)d9 



Act = TT 



TT-9 



As a second application of formula ( |1 . 2| ) we will consider the well studied example of 
QCD on a plane. We find that the structure of Wilson loops Wn{A) with large values of 
winding number n experiences qualitative changes as the area enclosed by the loop. A, 
passes through the critical value Api = 4. This phenomenon, which also can be viewed as 
a phase transition, is in fact a remnant of the Douglas-Kazakov transition occuring on a 
sphere of finite area A = n"^. 

We will use this result in section 3 to make some general comments about the impli- 
cations and meaning of the phase transition. 

The collective field theory approach to QCD2 is very powerful, at least on the cylinder. 
The Hopf action, as we shall discuss, propagates Wilson loops, or eigenvalue densities, along 
the cylinder. 

To discuss other geometries, in particular surfaces with handles, we have to know how 
to split the loops at a vertex. To this end we need the partition function for QCD2 on a 
"pair of trousers" -type manifold (fig. 2) (that is, a sphere with three holes). Moreover, it 
suffices to consider the limiting case when the area of this manifold vanishes. Indeed, we 
can always create the "pair of trousers" with a finite area by attaching cylinders to the 
zero-area "vertex". 

The properties of these string vertices shall be investigated in section 4. We find that, 
as — 00, the partition function of the QCD on a "vertex" vanishes in the zero-area 
limit, unless certain selection rule is satisfied. We find and prove this selection rule and 
discuss its implications for the string theory of two-dimensional QCD. 



See the discussion preceding (|l.lD . 
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Fig. 2: The vertex manifold with three boundaries, Ci, C2 and C3. 

2. Large QCD on a cylinder. 

The partition function of QCD on a cylinder (fig.l) can be evaluated exactly for any 



where A = g'^N is the large coupling constant and A the area of the cylinder. The 
summation is over all irreducible representations R of the gauge group, Xr(U) is the 
character of the matrix U in the representation R, and C2{R) is the quadratic Casimir 
operator of this representation. Since the answer depends only on the product of A and A, 
we shall set A = 1. 

If the gauge group is U (N), the representations R can be labelled by a set of integers 
+00 > > /2 > . . . > ^AT > -00. Then 



finite and equals |T^ [jT^ 




(2.1) 



R 



XRiU) 



detlle'^^'^*' 



(2.2) 



J(e^^«) 



where e*^*^; /c = 1, . . . , are the eigenvalues of the (unitary) matrix [/, and 



j<k 



Also, 
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In this section we will study the partition function (2J.) in the large N limit. This 
can be done in two ways. One way is to show that for large is dominated by a 

particular representation R. If this representation is found, it is possible to express all 
physical observables in terms of the indices {li} labelling R. The second approach is to 
construct an appropriate collective field theory describing QCD2. While the latter method 
is more physically transparent, practical calculations have so far been carried out using 
the first approach. For example, it is not at all obvious how to find the solution of the 
collective field theory which describes the strong coupling phase of QCD on a sphere. 

These two approaches are complementary. Below we will derive the equation ( |1.2| ) 
relating the dominant representation to the solution of the collective theory. This will 
allow us to provide a simple physical interpretation of the formulas for Wilson loops on 
a sphere, found recently by Boulatov, Daul and Kazakov and extend them to the general 
case of a cylinder. 

2.1. The dominant representation. 

Generally, infinitely many representations contribute to the partition function ( |2.1| ). 
However, as 00 the sum is dominated by a single representation. Indeed, for large 
the U{N) characters ( |2.2|) behave asymptotically as 

Xfl(C/)^e^'=[^-(^/^)'<^(^)] (2.3) 

with some finite functional S[py,(T]. In this formula it is implicit that we take the limit 
N ^ 00 assuming that the eigenvalue distribution of the unitary NxN matrix U converges 
to a smooth function! a'(^), 9 G [0,27r]. In addition, it is assumed that the distribution of 
parameters iji = h/N, which define the representation R, also converges to another smooth 
function py (y) , that we can call the Young tableau density. The functional S is, in general, 
not easy to calculate. However, in some important cases it can be found explicitly. 
As a result, 

Zn{Uc,.UcM) - J^exp Ar2|s[py(y),ai(^)] + S[py (y), a2(^)] 

R L 



^ The eigenvalues of a unitary matrix lie on the unit circle in the complex plane and can be 
parametrized as Xj = e*^^ . 
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where the bar denotes complex conjugation. In this formula the summation over all rep- 
resentations R can be thought of as "functional integration" over all possible distributions 
of the Young tableaux, pviy)- 

For large we expect that this "functional integral" is dominated by a saddle point, 
which can be determined from the equation 



d_ 
dy 



SE[py{y),a^{e)] ^5E[pY{y),a2{e) 



Spriy) Spviy) 
Then the free energy of the theory equals 

A 



(2.4) 



S[py(y),ai(^)] 



A 



+S[py(2/),a2(^)] -- j PY{y)y^dy--\ + O(iVO) 

where to simplify the notation we have used a shifted variable y = y — 1/2. It is also useful 
to keep in mind that the derivative of with respect to the area A (the "specific heat 
capacity" ) is given byS 

1 BFn _ 1 1 
~ ~24 

are always greater 



pYiy)y^dy. 



(2.5) 



dA 24 2 

Since the indices li are discrete integers, the differences k 
than one, and therefore the density pviy) must satisfy the constraint pviy) < 1- The 
saturation of this bound causes a phase transition for large QCD 0. 

To see how this phase transition occurs, let us consider the two-dimensional QCD 
on a sphere. We can regard the sphere as a particular case of a cylinder with boundary 
conditions Uc^ = = I (that is, (Ti{9) = a2{9) = d{9)). In this situation the large 
limit of characters is easy to calculate. Indeed, since 

XRil)=dR = l[ik-lj), 

we have 

S[py(y),5(^)] = ^ j pY{yi)pY{y2) ln\yi - y2\ dyi dy2 . 



Then (O) entails 



Py{u) du 
y-u 



Ay 
2 ■ 



The solution of this equation is 



py{y) = -^A-^. 



(2.6) 



^ Although pviy) by itself depends on the area, this dependence does not contribute to the 
specific heat due to the saddle point equation ( ^ 
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However, it obeys the constraint pviv) < 1 only if A < tt^. If A > then, instead, 
we must set pviy) = 1 in some interval y G [—6, h] and the solution becomes 



62 }y 

-nay ' \ y^' 

outside of this interval, for y G [6, a] . In this formula 



PY{y) = -^^/W^7W^^i(-^. -) (2.7) 
Tiav V a/ 



+1 

' ^ 2 7 l + xw2 ^(1-A;2t,2)(l_^2) ^ 



-1 



is the elliptic integral of the third kind and the constants b and a are fixed by 

k = b/a, k' = a/i — A;2^ 

a(2E(A;) - A;'^K(A;)) = 1, (2.9) 
aA^4K{k). 

Accordingly, the specific heat ( |2.5| ) is different for A < tc'^ and A > 7r2. In fact, it is 
possible to show that at the critical point A = tc'^ the third derivative of the free energy is 
discontinuous, indicating a phase transition. 

Obviously, a similar transition will occur on the cylinder. In this case the saddle point 
equation ( p.4| ) becomes very involved. To find its solution we will have to use an entirely 
different treatment of QCD2, based on the ideas of the collective field theory. 

2.2. The collective field theory. 

The large limit of two-dimensional QCD can also be studied with the aid of a 
collective field theory. The underlying idea of collective field methods is to make a change 
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of variables describing the theory so that the new variables (the "collective variables") 
would have a well defined large limit [ n | . 

In our case the partition function depends on the two sets of eigenvalues {aI^^ = 
e*^^ '10 = 1,..., A^} and {A^^^ = e*^^ ' |^ = 1, • • • , -^}- The specific eigenvalues do not tend 
to any limit as we take to infinity. However, the corresponding eigenvalue distributions 
ai(0) and o^iQ") do have a large limit. In fact, a\ and contain all the data one needs 
to evaluate the leading N ^ 00 asymptotics of Z^. 

Therefore, it must be possible to obtain a differential equation on Zn as a functional 
of a I and a2- To achieve this goal we will write down a differential equation satisfied by 
Zn as a function of the 2N discrete variables {^i^'*} and {0^^^}. Then we will change 
variables to cri(6') and a2{9)- 

As a first step, one represents ( |2.1|) in the form 



Zn(Uc,,UcM) 



1 



1 



E — " 

y,e{±^,±^,...} n [si 



det||e'^"^b''|| det||e-*^^^d" 



sm o sin " ^ " 



(2.10) 



with Vk = jj {Ik ■ 
equation 



where 



N-l 



) . It is now obvious that Z^ = Z^ e^'^^^ satisfies the differential 



1 



N 



C>2 



d ~ 



v{e^^^)z 



N 



n 



sm 



Ur 



(2.11) 



If Zn = expN^Fjsf, then Fn has a well defined large A^ limit F = lim Fn. The lim- 



N- 



iting functional F can be determined if we convert (|2.11| ) into an equation for Fn and 
then change variables from {^i^''} to o'i{9) in that equation. As a final result of these 
manipulations, it follows thall 



F[ai{e),a2{e)\A] = S[ai{e),a2{e)\A] 



sm 



dOdip-- I a2{0)a2{(p)ln 



e-if 
sm 



d6 d(fi 



(2.12) 



The procedure of changing variables, which is well known but nontrivial, is described in 
Appendix A. 
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where the functional is a solution of 



dS _ 1 
dA~2 



2tt 



d 5S V 7r2 



86 6ai{6) 



alio) 



(2.13) 



It is helpful to interpret (|2.13|) as the Hamilton-Jacobi equation for the classical Hamilto- 
nian 



271- 



H[a{e)Me)] = \ ja{e) 



(2.14) 



with A playing the role of time, the function (Ji{6) being the canonical coordinate and 



n(^) = SS/5ai{6) the conjugate momentum. The Hamiltonian (|2.14|) is called the Das- 
Jevicki Hamiltonian [|l4l . 

In fact, the required solution is easy to construct. To do this, we solve the Hamilton 
equations of motion 



da{e) 5H[a,U] 



dU{e) 5H[a,U] 



dt 5Ii{e) ' dt da{9) 

and pick the particular solution which satisfies the boundary conditions 

aie)\t^o = MO) 

cr {e)\t= A = a2ie) 



(2.15) 



(2.16) 



where, as before, cxi and a2 are the eigenvalue densities of matrices Uc\ and Uc2- Then it 
is possible to prove that S[ai, a2\A] equals the classical action calculated for this particular 
solutionli. 

The equations of motion which follow from ( p.l5|) are 



dv dv d f Tv'^ r 



dU 

de' 



(2.17) 



These are the Euler equations for a fluid with negative pressure P = — ^ . The solution 
we are looking for corresponds to the process where the density profile oiiO) evolves into 
the profile a2{0) during a time equal to A. 



If there are several such solutions, the one with the largest value of 5" must be chosen. 
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So far we have seen how this formahsm gives us a way of calculating free energy. 
However, in addition, the collective theory provides a natural and concise description of 
Wilson loops in the large N QCD. 

The application of the collective field theory to the study of Wilson loops is based on 
the following factorization property. Imagine we cut our cylinder along some contour C 
(fig. 1) into two pieces of areas Ai and A2 . We may ask how the partition function of the 
whole cylinder is expressed in terms of the partition functions for these two pieces. 

The answer to this question is obvious. Using (|2.1j ) along with the orthogonality 
relation for characters 

j dU xrAU)xr,{U^) = Sr,,r, (2.18) 

we easily see that 

Zn{Uc,,UcM) = j dUcZNiUc,,Uc\A,)ZN{Uc,UcM2)- (2.19) 

That is to say, to glue the two pieces into a single cylinder we have to set their boundary 
values, at the place we glue, equal to the same matrix Uc and then integrate over all 
such matrices. If we treat Z]s[{UcnUc2\A) as the probability amplitude for the process 
where Uc^ evolves into during the time A, then ( |2.19|) is just the convolution property 
satisfied by the Feynman path integrals. 

However, in the large limit the integral in ( p.l9| ) is dominated by a single saddle 
point. Thus the partition function reduces to the product of the partition functions of its 
two pieces. To see this, recall that the characters Xr{Uc) depend only on the eigenvalues 
of Uc, and so does Z^. Therefore, the integration in ( |2.19| ) reduces toi 

2iv({C^},{^f }|^) 



Substituting for Z^^ their expressions in terms of S using (|2.12| ), we obtain 



^N^S[<T^,a2\A] ^ j ^q{C) _ _ _ ^q{C) ^N^ {s[cTuoc\A^]+S[ac,02\A2]) ^ ^2.20) 

where (Tc{d) is the distribution of angles 0- over which the integration is performed. 
The integral in ( p. 20 ) is indeed dominated by a saddle point. We can even find this saddle 



^ The factor |^P(0'-^-')j is the jacobian which appears when we express the Haar measure dUc 
in terms of eigenvalues. 
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point exactly from the equation 
d 



(^S[ai, ac\Ai] + S[ac, 0-21^2! 



(2.21) 

d " 



(^S[ai, ac\Ai] + S[ac, ci2\M_ 



0. 



de 6aci0) 

Since S'[(Ti, acl^i] is the action along the classical trajectory connecting ai and o"2, the 
functional derivative of S equals the corresponding canonical momentum, 

d SS[ai,ac\Ai] 
de 6ac{9) 

where ^1(6*) is the velocity at the end of the trajectory, when a{t, 9) = ac(0). Similarly, 

d SS[ac,(J2\A2] 

with V2{0) being the velocity at the beginning of the trajectory connecting ac{0) to 02(0). 

Hence, the saddle point condition ( |2.21| ) states that vi{9) = V2{0). That is to say, the 
two trajectories ai ac and ac 0-2 can be joined so that the velocity is continuous. The 
resulting compound trajectory a*(t, 9) is just the solution of the collective field equations 
( p.lTp which describes the original cylinder of area A with the boundary conditions ai and 
a2- To find this trajectory, we must solve (|2.17| ) together with the boundary conditions 
( p.l6|) . Then the saddle point density ac{0) can be determined as a^{t = Ai,9). 

Using these arguments we can immediately evaluate any Wilson loop average on a 
cylinder. In general, a Wilson loop with the winding number n around the contour C 
equals 



W^{C)^(^TtUS 

(2 22) 

At large the integral in ( |2.22| ) is again dominated by a saddle point. Moreover, this is 
exactly the same saddle point which dominates the integral ( p.20|) . Indeed, the new term 
in the integrand, (l/A)Tr t/^ = (l/N) J2^=i 5 is finite as A ^ 00 and does not affect 
the position of the saddle point, which is determined by the balance of exponentially large, 
e-^ ^ , terms. As a result, we can calculate Wn{C) at large A evaluating (1/A) XljLi 
the saddle point, to obtain 

Wn{C) = J a,{t = Ai,e)e'''Ue. (2.23) 
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We have thus proved that the solutions of the coUective field equations have the physical 
meaning of yielding the dominant eigenvalue distributions which, upon Fourier transform, 
produce the Wilson loopJ^ for large QCD. 

2.3. The duality relation for QCD on a cylinder. 

Now we are in a position to find explicitly the representation which dominates the 
partition function. Let us first consider the symmetric case Uc-^ = Uc2- Later it will be 
obvious that our treatment easily gengeralizes to include all other boundary conditions. 
We will show that pviy) satisfies the equation 

TiPY{-Ticjo{e))=e (2.24) 



where cto(6') is the solution of collective equations (|2.17| ) taken at the middle of the cylinder, 
for t = A/2. 



To illustrate this relation let us reexamine the example of QCD on a sphere [|TT| . There 
cfiiO) = a2{0) = d{9) and the solution of the collective field problem ( |2.16|) , (|2.17|) is given 
by the self-similar evolution of a semicircular eigenvalue distribution 



a4t,e) = l^p{t)-^^-^ (2.25) 



where the scale n{t) changes in time according to 



Mt) = (2.26) 



Indeed, if we plug the ansatz ( |2.25| ) into ( p.l7|) we find that 



„2 

v{t, 6) = a{t)9, fi + 2ap = 0, a + + ^ = (2.27) 

whose solution is given by ( |2.26| ). At t = A/2 we have fx{t) = = 4/A and 

A \ 2 



ao{e) = a, (t = -, ^) = —y/A-9^. 



Using (|2.6|) we can check that, indeed, Trpy ( — 7Tao{9)) = 9 is satisfied. 



In general, the solution of collective field equations cannot be written down in an 
explicit form. However, these equations themselves can be integrated thus reducing the 



^ A different formula for Wilson loops on a sphere has been obtained by Boulatov, Daul and 



Kazakov [15|. Using the duality relation (|1.2| ) it is possible to prove that it is equivalent to ( |2.23 ). 
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problem to a single implicit equation on a^{t,6). To see this, let us introduce a new 
unknown function 

f{t,e) = v{t,e) + i'Kp{t,e). (2.28) 

Then the collective field equations ( |2.17| ) can be reduced to a single equation 

known as the Hopf (or Burgers) equation. Its solution can be found from the implicit 
formulall 

f{t,e) = u{e-tf{t,e)) (2.30) 

where the function fo{6) represents the initial data, fo{0) = f(t = 0, 9). 

We can use ( |2.30[ ) to relate a"o(6') to the density at the boundary, ai{9). Indeed, if 
the boundary densities are the same {(ri{0) = ct2(6')) then, by symmetry, the velocity at 
t — A/2 vanishes: 

Thus we can use 

fo{9) = f{t=^,9)=i7rao{9) 
as an initial condition in ( |2.3C1| ). After the time At = A/2 this /o should evolve into 

f{t = A,9)=vi(9) + i7iai{9), 

with a{t = A,9) = ai{9) and some vi{9) = v{t = A,9). Thus we obtain the equation 



ZTrao 



9-^{vi{9)+inai{9))] = Vi{9) + inai{9) (2.31) 



constraining ao{9) (and also vi{9)) as soon as <ti{9) is known. Although it would be very 
difficult to express ao through cxi in an explicit form, ( p.31|) contains just what we need to 
solve for the dominant Young tableau densiity pviy) in ([ 



To do it we will have to evaluate the functional derivatives of S[py, cti] in (2^). This 



can be done if we observe that the U {N) characters (^]5) can be represented as analytic 



continuations of the Itzykson-Zuber integral [|T6 



/.(^.B)./<if/e™<— '-^IJl^. (2.32) 



Since is a coordinate on a circle, f{t, 6) and fo{6) should be regarded as periodic functions 
of 6 in this formula . 
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In this formula A and B are arbitrary hermitian matrices, ak and bj are their eigenvalues 
and A(afc) = ni<j(^i ~ ^j) the Van der Monde determinant. Setting ak = bj = 9j 
and analytically continuing — > iak, we see that 

det ||e^"'=^^- II ^ J{e'^^) Xr{U). (2.33) 

Therefore, we can use the known expressions for the large A'^ limit of the Itzykson- 
Zuber integral [|T^ to find the functional S in ( p.3| ). In particular, if as ^ oo the 
distributions of {a^} and {bj} converge to smooth functions a(a) and then asymp- 
totically 



In{A, S) ^expA^2|^[a,/3] + ^ j a{a) a^da+]^ j 13(b) b"^ db 
~ 1^ j Oi{a) a{a) In \a — a \da da —- J (3 (b) (3 (b') In \b — b'\ db db' 



(2.34) 



with some functional 5" depending on the two distributions a and (3. 

Although there is an explicit expression for S[a, [3], we need to know only its functional 
derivative 

To determine it we construct a pair of functions \^G+{z),G-{z)^ such that 

G+{G.iz))=G.{G+iz))=z, 

Im G+ (a + iO) = na{a) , (2.36) 
lmG-{b + iO) = -np{b). 
Then 

U{a) = ReG+{a + iO) - a, (2.37) 

and, in addition, 

V{b) ^ -^S[a,(3] = -ReG_(6 + zO) +6. (2.38) 

If we want to apply these formulas to the calculation of characters, we must be able to 
perform an analytic continuation Ik ilk- This can be done by introducing an additional 
parameter t and rescaling Ik tlk- We will keep t in our formulas up to the very end of 
the calculations, where we will set t = i. 

In this case the density of points ykit) = = equals! 



Pt{y) = ]p{^). (2.39) 



^ Prom now on, we omit the subscript "y in pY{y)- 
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Therefore, we must distinguish between the derivatives with respect to p{y) and ptiv)- 
The quantity which enters the master field equation (|2.4| ) is the derivative with respect to 

p{y), 

(2.40) 



dy 5p{y) 

On the other hand, the constraints ( |2.36| ) will yield not U{y), but 

d 6 



Utiy) 



-S[pt{y),a^{e)]. 



dy Sptiy) 

However, these two derivatives are connected by a simple relation 



(2.41) 



U{y) = tUt{ty). 



(2.42) 



Indeed, if 



then 



Ht{y) 



Sptiy) 
5S[ptiy),ai{6)] _ 5 



S[pt{y),a^{e)] 



sp{y) 
s 



Sp{y) 



dz Spt{z) Ht{z) 



Spiy) 



dw5p{w)Ht{tw) = Ht{ty) 



where we made the change of variables w = z/t and utilized ( |2.39| ). Now it is easy to see 
that 

U{y) = ^Ht{ty)=tUt{ty). 



Therefore, using (|]3), (|]32D, (IjD, (|]3|) and taking into account that In J(e^^^ 
does not contribute to the derivative with respect to p{y)i we obtain 



dy 6p{y) 
= Uiy)+t^y 



t=i 



(2.43) 



tUtity) 



J t=i 



-y- 



Now we are able to produce the solution of (|2.4|) . If ai = a2, this equation takes the 
form0 

d dE[piy),a,{e)] _A^ 
dy Sp{y) 2 



(2.44) 



We will prove that if irpi^y) is the inverse function of 7rao(6') then (|2.44|) is indeed true. 



10 



We remember that y = y — \. 
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To this effect we exibit a pair of functions \^G^{z), G-{z)^ which satisfy the constraint 
for a{a) = pt(a) and (3{b) = ai{b): 



G+{z) = --Z + mpt{z), (2.45a) 
G_ {z) = -vi (9) - inai (9) (2.456) 



where vi is the velocity appearing in ( p.31|) . This would imply that, according to (|2.37|) 



Ut{y) = ReG+{y + iO)-y = --y - y 

and using ( |2.42[ ) we immediately see that ( p.44| ) is satisfied. 

The conditions (|2.36|) are easy to check. While the second and the third condition 
in ( |2.36D are satisfied by construction, the first one follows from ( |2.3lD . Indeed, we can 



rewrite (|2.31|) as 

TTao\9 + ^G-{9)] = iG-{9). 
Using 7rp(-7rao(6')) = 9 and pt{y)\t=i = jpi^)\t=i = -ip{-iy), we deduce 

9 + ^G-{9) = 7rp{-iG-i9)) = inpt{G-i9)). 

Thus 

-^G.{9)+i7rpt{G.{9))=9 

which, together with ( p.45a|) proves G+[G-{z)) = z, completing our construction. 

If the two boundary densities a\ and a"2 are not identical then the partition func- 
tion may or may not be dominated by a single representation. These two situations can 
be distinguished by looking at the collective field trajectory connecting a\ to ai- If at 
some moment to the velocity v{tQ,9^ vanishes simultaneously for all 9 G [0,27r], then the 
dominant representation exists and is determined byliil 

7rp(-7ra(to,6')) =9. 

On the other hand, if this is not the case the saddle point determined by ( p.4|) shifts to 
the complex domain and no single representation dominates Z\Uc^-,'Uc-2\-^- 



11 



Now to does not have to equal Ajl. 
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2.4- The Douglas-Kazakov phase transition on a cylinder. 

A phase transition similar to the Douglas-Kazakov phase transition on a sphere occurs 
in the large N QCD on a cylinder. 

Technically, this phase transition takes place when the range of values of the function 
Trpy (2/) fills the whole interval [0, 1]. Then in the strong coupling phase there are continuous 
intervals of y where pviy) = 1, while in the weak coupling phase everywhere pviy) < 1- 



Using the duality relation (|2.24| ) we can prove a more physical criterion which deter- 
mines the transition point. Indeed, since Trpyiy) and 7iao{0) are inverse function jii the 
range of values of npyiy) coincides with the domain of definition for nao{9). Therefore 
at the point of phase transition the domain of nao ranges from —n to n thus filling the 
whole circle. On the contrary, in the weak coupling phase this domain is restricted to some 
interval inside a circle and there is a gap where (Tq{9) vanishes. Thus this phase transition 
is precisely of the form found originally for the unitary matrix model [Q. 

Let us now look at the whole collective field trajectory a* (t, 9) . As we have seen from 
( p.lTp , a* evolves as the density of the fiuid with negative pressure. Since such fiuid tends 
to collapse the moment of time when its velocity vanishes corresponds to the maximum ex- 
pansion of the fluid. The density profile at this moment is, by definition, ao{9). Therefore, 
if the support of ao has a gap, this gap can only increase during further evolution, and the 
support of a*(t, 9) at any other t will have a gap as well. As a result, in the weak cou- 
pling phase the collective density exibits the gap at all times, while in the strong coupling 
phase there are time intervals when the gap disappears. This provides a criterion for the 
Douglas-Kazakov phase transition which is more general than the condition pviy) = 1- 
Indeed, it is applicable even when no single representation dominates the partition function 
and no pyiv) exists. 

Returning to the QCD on a sphere, we can now construct the collective field trajectory 
in the strong coupling phased. We are looking for a solution of the Hopf equation ( |2.2y| ) 
with the boundary conditions Im/(t = 0,6*) = Im/(t = A,9) = iTd{9). Surprisingly, the 
solution of this problem is quite complicated. Indeed, by the duality formula ( |2.24| ) the 
value of f{t = A/2, 9) = Tvao{9) is the inverse function of 7rpy(|/) which, in turn, is given 
by the elliptic integral (|2.7|). Since this solution describes the strong coupling phase the 
support of (Jo has no gap (fig. 4). 



The functional inversion in this context should be carefully understood since, strictly speak- 
ing, the inverse of npYiy) has several branches which should be chosen in an appropriate way. As 
a consequence, the domain of rrao is [— 7r,7r] and not [0,7r]. 

The trajectory in the weak coupling phase is given by ( p.25| ). 
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-a -b 



b a 



Fig. 4: The Young tableau density pviy) (left) and the collective field density 
o'o{9) (right) for the QCD on a sphere. In accordance with the duality relation the 
plots of TTpyiy) and 7Tao{6) can be obtained from each other by a reflection with 
respect to the axis x = y. 

However, the gap will inevitably appear as ao{9) evolves according to the Hopf equa- 
tion. After all, the final result of this evolution is the density cri(^) = S{9) with support 
consisting of a single point. 

Since our boundary conditions are even, so is o"*(t, 9) for any t. Therefore, the gap in 
the support of a* will be centered around 9 = n. In particular, we can find the moment of 
time when this gap just appears by looking at the solutions of cr*(t, 9 = tv) = 0. We will 
find that this equation has solutions only if t G [0, — to{A)] U [^A + to{A), A]. On the 
other hand, when t e]^A — to{A), ^A + to{A)[ the gap is absent and cr*(t, 9) ^ 0. 

Since the Wilson loops are nothing but the Fourier coefficients of 0"^<(t, 6*), we conclude 
that the structure of small loops (with area less than the "critical" value tcr = |^~'^o(^)) 
resembles the structure of loops in the weak coupling phase, even though we are in the 
strong coupling phase. On the other hand the loops with larger area are truly characteristic 
of the strong coupling phase. Physically, we can distinguish between these two types of 
loops by considering their behavior as the winding number n goes to infinity. Then at 
fixed area of the contour the strong coupling type of a loopliil decreases exponentially with 
n 



while the weak coupling type of a loop decreases at most algebraically, as 0{n p). The 
same is true for the general case of a cylinder. 

By definition, this is the loop for which the corresponding collective field density has no gap. 
Alternatively, the weak coupling type means that the gap is present. This should not be confused 
with the strong and weak coupling phases of the Douglas-Kazakov transition. 



Wr,{C) = C(e-"(^)"), a(C) > 



(2.46) 
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Continuing our example, let us find how the critical Wilson loop area depends on A. 
Obviously, to{A) vanishes at the Douglas-Kazakov transition, when A = tt^, and then 
grows as A increases. 

The evolution of a^{t,6 = n) is easy to study. We can use ( |2.3C1| ) with the initial 
condition fo{0) = inao {9) to obtainS 



However, due to the symmetry 9^—9 the velocity v{t, 9) vanishes at 6* = and 9 = tv. 
Thus f{t, tt) = z7ro"*(t, tt) and 

TTCTo 7T — i(^t — —^na^{t,n) =7va^{t,7i). 
Finally, using the duality (|2.24| ) we obtain an equation for z{t) — na^{t, tt), 

- i(t - ^y{t) = 7,pY{z{t)) (2.47) 

with pviy) given by (|3)- 

The solution z{t) we are looking for is a real function of t. Moreover, at all times 
Q < z{t) < h. Indeed, dX t = Iq = A/2 the expansion of our hypothethic fluid is maximal 
and so is its density at 6* = tt. When the fluid collapses towards 6* = 0, a(t, tt) decreases to 
zero. 

Mathematically, if z G [0, h] 

npyiz) = i — y^{a^ - z^){h^ - z^) Hi f-^, k). (2.48) 
az \ z'^ J 

Since in this case x = — ^ < —1, the elliptic integral ni(x, /c) develops an imaginary part. 
This happens because the poles of the integrand at = — ^ are now located within the 
integration region. Using the rule ^ = V\ — Z7r5(a;), we get 

+1 

ni(x, k) = H , ^ - (2.49) 

2_(l + ^«' ^{\-k^u-){\-u^) y/|x|(l + f)(l + i) 



Substituting this into ( |2.48| ) and using k = h/a, we see that 

+1 ^ 

npY(z) = TT H a/ (a^ — z'^)(b'^ — z"^) + — 

az^ ' 7 (1 



(1 + xu^),/{l-k^u^){l-u^) 



Since we set the initial condition at to = A/2, the time of evolution in this equation is t — | A, 
not t. 
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Therefore ( |2.47| ) reduces to a single real equation 



+1 



(t- -)z = - — J(a^ - z^W - z^) i r2 > (2-50) 



-1 



It is easy to see that 2; = is always a solution of this equation. However, for certain 
t it has another solution which can be determined from 



t-^ = --^/{a^-z^){b^-z^)-f 
2 a J 



+1 

du 



(z2 - U^b^)^{l-k^U^){l 



Such z{t) decreases as t — ^ deviates from zero, so that when t goes outside of the interval 
]\A — ro(A), + To(A)[ this solution ceases to exist. Therefore, at t = + tq{A) we 



have z{t) =0 and 



To (A) = -6 lim 

2— >o 



+1 

du 



-1 



Using the x — > asymptotics 



+1 1 



/du , 9 /" u du Ox / X 
= 2k^x / , + Oix^) 2.51 
^ (w-x)v/(l-A;2w2)(i_„2) y v^(l-A;2w2)(i_^2) 



we conclude that 



^.^y _ 2A;2 } u^du 



and finally 



b 

1 /" A2 r/A 2 

-°'^'^^/ v(o--A^)(>--A^) ^><'^W-EW). (2.52, 
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Fig. 5: The critical Wilson loop area tcr as a function of the size of the sphere. 
As A ^ oo (this corresponds to the QCD on a plane) tcr converges to a finite limit 



The dependence of the critical Wilson loop size tcr = \^ ~ ''"o(^) on ^ is shown in 
fig. 5. Most remarkably, tcr tends to a finite limit equal to tcr = 4 as A goes to infinity. 
Thus even in the QCD on a plane the eigenvalue distributions for small contours have a 
gap which disappears as the area enclosed by the contour passes through 4 That is 



to say, the master field for QCD on a plane contains some trace of the Douglas-Kazakov 
phase transition and of the weak coupling phase. 

Since the asymptotic property tcr 4 will be important to us later on, let us prove 
it. Making the substitution A = hcosip in (|2.52D we get 



ro{A) 



■k/2 

2 f cos^ (f dip 



a 







+ sin If 



where = ^ — 1 = At A — > oo, as a consequence of (^), k''^ ^ 16e~^/^ and 
thus = 16e~"^/^ + . . .. In this limit we can break up our integral into two parts and 
construct the asymptotic expansions for both of them: 

7I-/2 e 7r/2 

cos^ (fdip f dip ^ r cos^ (p ^ 



^ ^ ^ e 

= arsinh- + Incotan - - 1 + 0(e) + OU^) 
K 2 

= _1 + In 1 + . . . = -1 + 4 + 0(e-^/4) 
n 8 



21 



In these formulas we chose an auxihary smaU quantity e so that k ^ e ^ 1. Then at the 
end e drops out of the final result, leaving us with the desired asymptotics. 
Now from (^) a = i + 0{e-^/^) and thus finally 

ro(A) = |-4 + 0(e-^/4) 

leading to 

ter(^) = 4 + 0(e-^/^) (2.53) 

as claimed. 

Let us now demonstrate that the asymptotic behaviour of Wilson loops for large 
winding numbers n is related to the Douglas-Kazakov transition. The most convenient 
way to do this is to continue considering the planar case. For one, as we saw, the QCD on a 
plane retains some trace of the weak coupling phase (the gap in the eigenvalue distribution 
for loops with area less than 4). On the other hand, all Wilson loops on a plane can be 
evaluated explicitly, making the analysis very simple. 

Indeed, the Wilson loops on a plane equal ^ Jl^ 



Wr^iA) = {^TrUS) = V^4^2i(An) (2.54) 
where L^j^}_i{x) are the so-called associated Laguerre polynomials 

n — 1 

Using this expression we can investigate the behaviour of Wn{A) as n — oo at fixed A. 
We will see that VF„(A) decays exponentially in n if A > 4, but as a power of n if A < 4. 
To prove this we need to know the asymptotic behaviour of the Laguerre poynomials as 
X = An oo. This can be found by treating the integral in (|2.55|) using the saddle point 



methods. The asymptotics of this integral is governed by the minimum of the functional 

$(t) = -t-lnfl + i). (2.56) 

This minimum is reached at 



n \ t 



2V4a; 2V4A 



This calculation is remarkably similar to the calculation of instanton contributions in the weak 
coupling phase which are expressed as Laguerre polynomials as well Q . This is not suprising since 
these are different manifestations of the same phase transition. 
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If A > 4 then both of and t- are real. Inn this case it is possible to show that the 
dominant contribution comes from the region around t = and equals 

L«,(x) ^ ^ ^ e--^^^^) = -^(l- (2.57) 



'27m|$"(t+j| 
where we have introduced a special function 



7M^v^-^ln | + ^ = 2v^^il-^>0. (2.58) 
Thus at A > 4 

so that as n — > oo, VF^(A) decays exponentially with the index a{A) = ^7(4/A) > 0. 

On the contrary, if A < 4 the saddle points t+ and t- are complex conjugate and 
make comparable contributions. Then 

4-i(2^) - (-)""' ? 2Re {e-"^(*+)-^"/n 

/27rn|$"(t+)| ^ ^ 

where r(x) = a/x — 1 — xarctan-v/a^ — 1- As a result, for A < 4 the Wilson loops 



(2.60) 




-"(-'-&(l-0"'-{xKl)n} <-^' 

do not decay exponentially. Rather, they oscillate with n and the amplitude of oscillation 
decays as n~^l'^ . Such behaviour is similar to what we encounter in the weak coupling 
phase on a sphere before the Douglas-Kazakov transition occur 

To summarize, the phase transition in continuum two-dimensional QCD exibits itself 
through physical observables in several ways. First, the density of eigenvalues for Wilson 
loop matrices develops a gap in the weak coupling phase. Second, the asymptotic behaviour 
of Wilson loops with large winding numbers is different in the two phases. Finally, the 
Young tableau of the dominant representation also can be expressed in terms of physical 
observables using the duality formula ( |2.24| ). 



Let us emphasize that the absence of a gap and the exponential decay of Wilson loops are 
not automatic consequences of each other. In fact, the n ^ oo behaviour of IVn(A) probes the 
degree of smoothness of the eigenvalue distribution. Since when there is no gap the distribution 
is infinitely smooth the Wilson loops decay faster than any power of n. On the other hand, if 
we do have a gap, then at the edge of the support of the eigenvalues higher derivatives become 
discontinuous, and the loops decay powerlike. 



23 



3. Instanton contributions to the Wilson loops. 

As we saw in the previous section, the Douglas-Kazalcov phase transition is associated 
with the development of a gap in the eigenvalue distributions for Wilson matrices. In this 
section we will argue that this phenomenon, very much like the original Douglas-Kazakov 
transition on a sphere at Acr = tt^, is induced by instantons. 

Indeed, it is known that the partition function of QCD2 can be repesented exactly as 
a sum of contributions localized at the instantons, that is, the classical solutions of the 
theory As it turns out, in the weak coupling phase, at Acr < tt^, the dominant term 



comes from the expansion around the perturbative vacuum, ^^(x) = 0, while in the strong 
coupling phase, Acr > tt^, a certain nontrivial instanton configuration will dominate. 

A similar analysis can be carried out for Wilson loops. We can anticipate that the 
average of a Wilson loop will be given just by its classical value as determined by the 
dominant instanton configurationlii modified by quantum corrections. The easiest way to 
determine these is to perform a Poisson resummation in the original formula for Wilson 
loops in terms of the representations of the gauge group. On the sphere, for a simple 
contour which divides it into two patches of areas Ai and A2, the Wilson loop is given by 



X 



1 

N 



(3.1) 



In this formula R and S label all irreducible representations of the gauge group, dn and 
ds are their dimensions and x{U) are the characters given by (|2.2|) . (|3.1|) can be derived 
easily by gluing two discs of areas Ai and A2 along the loop. 

Labelling R and S by two sets of integers, > 4^'' > . . . > and > 4^"* > 
. . . > l^j^^ and transforming J dU to an integral over the eigenvalues of U we get 



WniAi, A2)Z{A^ + A2) = A(4^^) A(/f ))e-^^^ 



'1 ^'2 ^--^N 

'1 ^'2 ^■■■^'iv 



(3.2) 



X 



r / 1 \ 

/ d6i ...d6N( j^Y^""'^'' ) det|| e^'^'^^^^ll det|| e" 
^ k=i ^ 



By instantons we mean all possible solutions of the classical Yang-Mills equations even if 
they are unstable. 

This is the value we obtain if we substitute the classical Af_i{x) for the dominant instanton 
into trU = tr Vexp ^ (x) dx^ . 
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with A{xi) = Y[i<ji^'i ~ ^j)- Now we can expand the determinants and use the antisym- 
metry of A(/) with respect to the permutation of Zj's to remove the ordering restriction on 
the Vs. We obtain 

1 ^ 

WUA^,A,)Z{A^ + A,) = -Yl E A(/fV(^f ) 

k = l ;(S) 

i ' 3 

unrestricted 



„ N N 

p=l q = l 



1 ^ 



fe=l ^(H) 



where now the summation is only over I - , the Ij being determined by /) — +^'^j,fc- 
That is to say, l^^^ = Ij^^ for aU i except for i = k, when Z^'^-' = l^^^ + n. 

To represent (|3.3|) as a sum of instanton contributions let us transform it using the 
Poisson resummation formula, 

+ 00 +00 

^ /(ni,...,nAr) = ^ F(27rmi,...,27rmAr) (3.4) 

rji ,...,njv = — tx) mi ,. . . ,mjv = — CX3 

where 



+ 00 

F(pi,...,p^)= y /(xi,...,XAr)e^(P^-i+-+P-^-)dxi...dx^. 

— oo 

First we recall that, up to an irrelevant uniform shift of l^s by (A^ — l)/2, 

TV 
i=l 

Thus 

and we can represent our Wilson loops in the form 

Wr.{A^,A2)Z{A^+A2) 



1 ^ 



N 



k = l 
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where we have introduced a new function of variables {Ij} 



<p{{lj}) =A{lj)e 4iv2.p=i'p 



(3.5) 



A = Ai + A2 being the total area of the sphere. 

To apply the Poisson formula ( |3.4[ ) we have to find the Fourier transform of the 
function under summation: 

F(pi,...,pAr) 

+00 

= j dxl...dxNV^{{x^})ip{{x, + 5,,kn})e-^^'^'-'^'^'"' e^(pi-i+-+P^-iv) 

—00 
+00 



J dxi... dxN <^({x J) ip{{x, + 5,,,n}) e^(pi-i+-+P— ) 



where pj = pj + {in/2N)(A2 — Ai)di^k- Recalling that the Fourier transform of a product 
is a convolution of the individual Fourier transforms we get 



+00 



F{p^,...,p^)= I dy^...dyN^({^^^])i^(^^^J^\]e-^^^^-y-^-/' 



(3.6) 



where is the Fourier transform of ip[{x}). 



+00 

/N 2 
rfxi...rfa;7v¥?({x,})e*(Pi"i+-+^^"^) = C7vA(pi) e"^ ^-1 (3.7) 



with some A^-dependent constant which will cancel later onE2l. The phase factor in 
is due to the shift of Xk ^ Xk + n in (/?({xi + di^k'n}). Thus 



F{pi, . . .,pn) = CNe~ 



/9 -2 



+00 



dyi. .. dyN 



■i<j 



1^ _jvy^^ „,2 



To derive this one represents A(Zj) in (3^) as a Van der Monde determinant and performs 
the integrations exphcitly. The answer can be simphfied to give ( |3.7| ) if we remember that one 
can add rows of a determinant without changing it. See H for the details. 
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where pij = pi —pj, Uij = Vi — yj and Cn is another constant, Cn = C'^/2^^^~'^\ 

Combining the pieces we finaUy obtain the foUowing representation for the Wilson 
loops: 

W^{AuA2)Z{Ai + A2) 

N +00 



^ N +00 



Er=i "'j p-27rmmfe^2M f.n'' (Ai-Asf /8AN (3- 



k=l mi,...,mjv = — 00 

N +00 , 

_2tt±I 

g A ^j = l J g g. 

k=l mi,...,mjv = — 00 

+00 

/N V^-^V 2 



where = rfn — rhj and mj = mj + in5j^k{A2 — Ai)/4'kN. 

This formula has an interpretation in terms of instantons. Indeed, the instantons of 
QCD2 can be labelled by integers mi, . . . .rriN and have the action 



6mst(mi, . . . , mN} = } 
The corresponding field configuration is just a collection of Dirac monopoles 



/mi^J(x) 

m2^°(x) 





V 

where = ^^(O, ((>) is the Dirac monopole potential. 



\ 







(3.9) 



(3.10) 



^^(e,(^) = o, ^«(e. 



1 — cos G 



O and (/) being the polar (spherical) coordinates on 5"^. We can see that the terms 
^■KznmkA^IA (|3]8|) represent the classical contribution of the field configurations ( |3.10| ) 
to the Wilson loop while the term 



+00 

-00 «<J 



2^ Z^j = l 



(3.11) 
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represents quantum corrections due to the fluctuations around the instanton. Thus our 
final result can be represented in the form 



1 r 1 ^ 

^ ' all 

instantons 



(k) 

N ' " 

k=i 



In fact, if we consider (|3.8|) with n = 0, we will reproduce the instanton representation 
for Z{A) constructed by Minahan and Polychronakos 

+00 

Z= J2 «^(mi,...,m7v)e-^^"^i'-'""^) (3.13) 

mi,...,mjv = — 00 

where 

+^ TV AT ^ ^iv 2 

w(mi,...,mAr) = / J] 11 (47r V, - yj) e ^ ^» (3.14) 

is the statistical weight of the configuration (mi, ...,111^). 

In the large limit we expect that the sum in (|3.13[ ) is dominated by a certain 
configuration (m^'^^ , • • ■ , ^tv^) depending on the area A. (For example, in the weak coupling 
phase the dominant configuration is m^i^ = . . . = m^"* =0.) As 00 for any such 

configuration the integrals (|3.11|) and (|3.14| ) are dominated by some saddle points in y^'s. 
Moreover, the saddle points for these two integrals are exactly the same, as they are 
determined by the balance of e'-^*^-^^) terms which are common for both integrals. Thus, in 
the large limit, we can determine ci!'^ /w by substituting the values i/i = yf^ into the 
ratio of the integrands: 



Note that, if Ai = A2 = A/2 (that is, if we are considering the Wilson loop around the 
equator) then rrii = fhi and 

Ak) 
W 

Therefore in the leading large order 



N Jk) . N 



^in{-iTm^.+y^°^ /2) 



fc=l k=l 

Thus we see that the large A^ density r{r]) describing the distribution of numbers 77^ = 
y^^^ /2 + Tinik is nothing but the eigenvalue distribution for the equatorial Wilson loop, 
o"o(^). In the weak coupling phase, where all = 0, this identity is easy to check directly. 
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Let us also mention that the distribution of {'mf'\ . . . , m^'') for the dominant instanton 
configuration in the strong couphng phase can be determined from the formula 

+00 +^ +^ 

J2 Pme--'^'' = fdh expj-^ [/ PyM^ ^Ah]\ = l-2b+ I dhe^'^^-^'^^ (3.15) 

— —0 

where q G [— tt, tt], is the probability to find an integer m among the set {w!{ \ . . . , rrS^} 
and pviy) is the Young tableau density given by (|2.7|) . This can be derived using the 
Poisson resummation in a way similar to ( p.8|) . 

As we can see from (|3.12|) , the quantum corrections make important quantitative con- 



tributions to Wilson loops. However, the physics of the phase transition can be understood 
even if we neglect them. Consider the classical contribution to the Wilson loop, 

1 ^ 

W''\Ai, A2) = —Y1 e-2™""^^^/^. (3.16) 
k=i 

The effective eigenvalue distribution for such a loop will have a gap if the angles 9^ = 
2'KmkA2/ A are restricted to a certain domain inside of a circle, \9k\ < ^max < for all k. 
Since 6'niax = 27rmjnax^2M5 we get the condition on the area of our loop 

Ao < 



2?Tlmax 



As A2 increases this inequality may be violated and the eigenvalue gap will disappeaiS. 
The exact value of A2 when this will occur will be smaller than A/2minax since quantum 
fluctuations cause some additional widening of the eigenvalue distribution. 

In general, for A2 = we always have 6k = 0, even in four-dimensional QCD, since 
the corresponding Wilson matrix is the identity. On the other hand, as the area of the loop 
becomes very large the gauge fields at the distant points of the contour are uncorrelated 
and the distribution of angles becomes uniform: cr{6) = l/2n. This is also true in the 
four-dimensional case. We can see that while zero-area loops have an eigenvalue gap, the 
infinitely large ones do not. This might be a piece of evidence for a large phase transition 
in QCD in any dimension. However, we also saw that the quantum corrections tend to 
broaden the eigenvalue distribution. In fact, it might happen that for any small nonzero 
area the distributions develop infinitesimal tails and do not have a gap anymore. The 



On the other hand, in the weak coupling phase ah mk = and the inequahty A2mmax < A/2 
is always satisfied, consistent with the fact that there is always a gap. 
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nontrivial fact about QCD2 is that such a phenomenon does not occur there. The reason 
is that the Hopf equation, governing the evolution of the eigenvalue densities, essentially 
prohibits the development of such smooth tails. Rather, its generic solution has a (moving) 
edge and vanishes as a{9, t) ~ \/0'^-^^^At) iisar it. Whether or not the situation is similar 
in higher dimensions is unclear. 

One could contemplate analysing the 4d Wilson loops for small area by perturbation 
theory, since in this case the asymptotic freedom applies. However, to obtain the eigenvalue 
distribution one would have to evaluate the Wilson loops with arbitrarily high winding 
numbers n. This corresponds to multiplying the charge flowing around the loop by n, and 
for high n this takes us out of the domain where we can rely on perturbation theory. 

Let us emphasize that the phase transition in the eigenvalue densities of the Wilson 
loops does not yet imply that there should be a singularity in the free energy of the 
theory. For example, this is the case with QCD on a two-dimensional plane. But such a 
phenomenon might imply a phase transition for the free energy of QCD in a box of a finite 
sizeS (say, for a sphere in the case of QCD2). In general, such transitions would mean 
that the string representation of QCD4, if constructed, will have only a certain domain of 
validity. On the other hand, it may offer an explanation how the short distance behavoiur 
of QCD which is described in terms of particles, agrees with the string description. 



4. Exact critical area for the phase transition on a disc. 

As another application of the duality formula ( |2.24| ) we will calculate the critical area 
for QCD on a disc. A disc can be regarded as a particular case of a cylinder with one of 
the boundary holonomy matrices set to identity, = that is to say ^2(6*) = 5(6*). 

Quite remarkably, in this case the boundary value problem for the Euler equations 
( P.16|) , (|2.17|) can be reduced to a Cauchy problem and then solved explicitly using ( |2.30|) . 



To this end we need to know how to specify the velocity v{t = 0, 9) = vi{9) at the 
beginning of the trajectory so that at the end, when t = A, the second boundary condition 
(|T|) is satisfied, a{9)\t=A = (^2(9) = S{9). 

While such a problem cannot be solved explicitly for an arbitrary (72{9), in the par- 
ticular case a2{9) = d{9) the necessary velocity is given by 



22 



See also the recent paper by Douglas [21| who draws similar conclusions. 
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To see why this is true let us recaU that the boundary value problem identical to ( |2.16|) , 
( p.lTp occurs in the calculation of the Itzykson-Zuber integral. The correspondence be- 
tween these two problems is established by cri(^) ^ 0(6'), 0-2(0) ^ (3{9) (see ( f2.36D and 
the analysis in the Appendix B). In the case 02(0) = 6{9) we would have /3(6) = 6{b) which 
means that the argument B of the Itzykson-Zuber integral is simply zeroil. Then the 
integral itself is obviously equal to one, In{A,B = 0) = 1. On the other hand, we can 
apply the asymptotic formula ( p.34|) to obtain 

S[a{a), (3{b) = 5{b)~\ = —- J a{a)a^da+- J a{a) a{a) ln\a — a'\ da da' . (4.2) 

The functional 5" is the action along the classical trajectory connecting a{a) and /3(6) 
during the unit time interval, t = 1. The velocity at the beginning of the trajectory can 
be found as 

, , d SS\a,f3] f a(a')da' 

^ " = TT X f V = -« + f , ■ 

oa da[a) J a — a' 

If the time interval is not equal to t = 1, this formula can be easily generalizedii 



, , a f a(a') da' 
v[a) = -- + 



t J a — a' 
finally proving ([4.1| ). 

Since the Euler equations ( p.l7| ) are equivalent to the Hopf equation ( |2.29|) , we can 
now set up the Cauchy problem for ( |2.29| ) with the initial condition 

n^ . ^ fai(e')de' e raAe')de' ,^ , 

f{t = 0,9)^ fo{9) = -j+j- + ^""^^^^^ ^ " A + y e-e' ■ ^^-^^ 

Then the solution is determined by ( |2.30|) which for this particular case translates into 

This equation can be used to determine fit, 6) and, consequently, a{t,6), for any specific 
<7i . However, the problem of the Douglas-Kazakov phase transition can be solved without 
doing so. Indeed, as we proved in the previous subsection, in order to test for the transition 



Note the important distinction from the matrix Uc2 1 which is the identity matrix for (J2{0) = 

5{e). 

To do this one considers the integral lN{A,B\t) = J dU e"r'^(^'^^'-^ \ Its asymptotics is 
obtained by replacing ^ J a{a) da ^ J a{a) da and | J (3{h)h'^ dh ^ ^ J /3(6) db in the 
formula (|]34|). 
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it is sufficient to find out whetlier or not tlie support of a(t, 9) develops a gap at any time 
t G [0, A]. Tliis allows one to obtain an explicit formula for tlie area of the disk when the 
transition takes place. 

To simplify the analysis, let us consider the case of a symmetric ai (9) , so that ai {9) = 
ai {—9). If the support of ai is the whole circle = [— tt, tv] then in the process of evolution 
this support must always develop a gap. Indeed, the final result of such evolution is a delta 
function with support which covers only a single point ^ = 0. 

Such picture is characteristic of the strong coupling phase. Therefore, in this situation 
we never have any weak coupling phase and the transition is impossible. A different picture 
arises if the support of o"i covers only a part of a circle [—6, b] with b < n. Then the Euler 
evolution can cause this support to expand with t before it starts contracting to ^ = 0. If 
this expansion stops before the edge of the support, reaches tt then the system stays 
in the weak coupling phase. Otherwise the transition to the strong coupling phase will 
occur. 

At the moment tc when b(t) reaches its maximum /(tci 9 = b{tc)) = 0. In- 
deed, (7"(tc, 6')|e=6(tc) = because b{tc) is the endpoint of the support of a. Moreover, 
v{tc, b{tc)) = because exactly at t = tc the support has just stopped expanding and the 
velocity at its edge is zeroiH. Thus we may determine 9c = b{tc) using ( |4.4| ) with / = 0: 

9c fai{9')d9' 



A J 9c -9 



The critical A = Acr corresponds to = giving 

ai{9) d9 



Act = TT 



TT 



9 



(4.5) 



Note that the integral in ( |4.5| ) is not a principal value integral. Rather, it is an ordinary 
integral which has no singularities because 9 = tt lies outside of the support of ai. 

Like the duality relation formula (^]^) can be checked in a number of exactly solvable 
cases. For example, QCD on a sphere corresponds to (yi{9) = 6{9) so that / "^^^J*/^ = ^■ 
Thus yields Acr = tt^ reproducing the original result of Douglas and Kazakov. 

The same technique can be used to obtain the transition area for a cylinder whenever 
we can evaluate the large limit of the corresponding Itzykson-Zuber integral ( p.34|) with 



the eigenvalue distributions a = ai and (3 = a2- This way we can generalize ( |4.5|) to the 
case of a flat a"2, that is, 

a2{9) = { 2c2 

otherwise 
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Still, at t = tc the velocities inside of the support do not have to vanish. 
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Then the C2 — ^ Hmit would correspond to a2{9) — ^ '^(^)- 

('2) 

For such a2 the eigenvalues 6j are distrubuted uniformly over the interval [— C2, C2] 



and the determinant in the Itzykson-Zuber formula (|2.32[) reduces to a Van der Monde 
determinants 



det 



^2c2ak{j-^) _ g- 



^^2 2^fe=i«Met (e^^^afe^J =e-^^2Z^fe=i"^ A(e2^2afej 



exp 



TV ^ TV 



A;=l 



Thus the large limit is easy to evaluate with the result that the action functional 5" 
equalj^^ 

s = -^j ai{6) e^d6-j j ai{e)ede 




2co n 2c9 /i/ 1 

-^^-e^^ dode'. 



(4.6) 



Therefore the initial Hopf velocity solving the boundary problem ( p.l6|) can be written 
down as 



v^{9) 



d 5S 



^C2 



d6 6ai{6) t 
ai{e') dO' 



9 + C2 ^2c2 f ai{9')d6' 



t J i_Q^{0'-e) 



(4.7) 



t t j tanh[f (^-^')] 
so that the initial function /o in (|2.3UD is 



ai{9')d9' 



t t j tanh[f (6'-6'')] ' 
As before, the critical area Acr can be found from 



/o(^c = 7r)|^^^^ = 



giving finally the formula 

ai{e) de _ Ti_ 

tanh[^(7r-0)] ~ V2 

We keep in mind the correspondence ^[^^ ak, O^p hj and a{B) 
See footnote 24 for the explanation of t-dependence. 



(4.8) 
(4.9) 

ai(0), ^{9)^a2{e). 
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In the limit c — this equation reproduces ( [4 .51) . On the other hand if ai also is a flat 
distribution, (Ji{6) = ^ for 16*1 < ci then gives us 

tanh-^^^ tanh— ^ = tanh^^i^. (4-10) 



5. Two-dimensional QCD on a vertex manifold. 

Another nontrivial problem is the large QCD partition function on a vertex manifold 
(flg. 2), or pants diagram. In a sense, the cylinder and the "vertex" are the only nontrivial 
manifolds we have to consider. Any other two-dimensional surface can be constructed 
by gluing together an appropriate number of vertices and cylinders. Then, according to 
( ^.19| )-( p.22D the partition function on such composite manifold is just the product of 



partition functions for its constituents. Moreover, we can restrict ourselves to the limit 
when the area of the vertex manifold goes to zero. Indeed, we can create a vertex of any 
finite area just by attaching cylinders to a vertex of infinitely small area. In this case the 
partition function of QCD on a vertex becomes a certain functional of the three eigenvalue 
distributions for the boundary matrices^ Uc^^, and Uc^- 

This functional should be somewhat similar to a delta-function. For one, in the case 
of a cylinder if A — 

Zn{UuU2\A) ^ 5^Xi?(f/i)Xi?(t^|) =4i(t^i,t^2) (5.1) 

R 

where 5c\ is the delta-function on the set of conjugacy classes of the gauge group. That is 
to say, 5c\ equates the sets of eigenvalues oi Ui and U2- We can represent it in terms of 
the usual delta-function on a group manifold 8{U): 

5c\{Ui,U2)= j dU6{UiUUlu\l) (5.2) 

with / the identity matrix. 

It is useful to think of this function as analogous to the vertices of conventional field 
theory. There, with every vertex there is associated a delta-function of the total incoming 
momentum. Obviously, something similar should appear in our case as well. We expect 
that the structure of the partition function of QCD on the vertex manifold will have two 
ingredients. One is a selection condition (similar to the condition that the sum of all 
incoming momenta is zero) and the other ingredient is the actual value of the partition 



Throughout this chapter these matrices will be denoted simply as t/i, U2 and U^. 



34 



function when this condition is satisfied (the counterpart of the vertex coefficient in field 
theory). 

To find the selection condition we start with the exact expression for the partition 
function of QCD2 on the vertex manifold of a finite area A 

2iv(C/i, Us\A) = V XRiUi)xRiU2)XRiUs) ^-^c.^n)_ (5 3) 

R 



As A ^ 0, 

^ \^ Xr{Ui)xr{U2)xr{U3) _ ^ jj jj 



R ^« 



Using the formulas 



j dVxR{AVBV^) = j-XR{A)xR{B) 



(5.4) 



R 

we can rewrite this in the form analogous to ( |5.1|) : 



Zn{Ui,U2,Us)=J2 I dViXRiViUiV^U2)xR{Us) 
dVi dV2 5{ViUi V2 U2 Us, I). 



(5.5) 



Since the characters in ( p.3| ) are class functions on the group we can assume that C/i, U2 
and Us are diagonal matrices. 

To evaluate the integral in (^]^) we need to find matrices Vi and V2 that satisfy the 
equation V1U1V1V2U2V2US = I and then compute the Jacobians necessary to integrate 
out the delta-function. 

First, we rewrite ( |5.5| ) as 

Zn{Ui,U2,Us) = I dV^dV25{V^{U^{VlV2)U2{VlV2)^)vlul) 

^ , (5.6) 
dVx dM5 {Vi {UiMU2M^)V^, . 



Now, we can perform the integration with respect to Vi, keeping M fixed. The delta- 
function will pick up the specific value of Vi = Vi^^ such that vj^^ diagonalizes the matrix 
Q{M) = UiMU2M^ . Thus, after the integration over Vi the delta-function in ( |5.6D reduces 
to a delta-function equating the eigenvalues of Q{M) and U^. Indeed, if v}'^^ Q{M)Vi^'^'^ = 
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A where A is diagonal then in the smaU vicinity of V-[^^ we have Vi = (I + vi)v}^^ with 
some antihermitian vi and 

5iViQ{M)vl, Ul) =5{A+ [vi, A], C/j) 

AT 

= llS{Re[vi,AU)S{lin[vi,AU) llS{Ap ~ {ul)p) 

p<q p=l 

= n ''''°'7'"'1';""''°' n ^(A. - (ci).) 

p<5 ' ^ ' p=l 

where Ap and (f/3)p are the eigenvalues of the (diagonal) matrices A and U3, respectively. 
Thus, up to a A-dependent coefficient c(A), 

N 

Zn{Ui,U2,Us) ^ dMc{A{M))ll5{Ap(M)-{Us);) (5.7) 

p=i 

where Ap(M) are the eigenvalues of Q{M) = UiMU2M'^ . To estimate this integral we 
should find M = Mq such that Q{Mo) has the same eigenvalues as t/J and compute 
the Jacobian J{M) = det|| ^g^^^ \\^ ^_]^- The partition function will be inversely 
proportional to this Jacobian. 

If this Jacobian vanishes the partition function blows up. Intuitively, in this case 
the volume of integration over dM which contributes to (^]3) much larger than in the 
generic case. The most singular situation is when not only the Jacobian J(M) bot the 
whole matrix ^gl^^^ vanishes. Then any small shift away from M = Mq does not violate 
the condition imposed by the delta-function in (|5.7| ) and the integration volume will be 
the largest. The condition on Ui, U2 and C/3 under which this occurs will be an analogue 
of the condition Yll=i Pi = ^ imposed by the vertices in field theory. 

To find this condition we must calculate the variation of eigenvalues 5Ap(M) — 
Ap(M + dM) — Ap(M) to the first order in 5M and require that such variation vanishes. 

But Ap(M) are the eigenvalues of the matrix Q{M) with the variation 

5Q{M) = Ui5MU2MI - U1M0U2MI5MMI. 

If 

v/°^Q(Mo)y/°^^ = y/°ViMof/2MoV/°^^ = A(Mo) (5.8) 

is diagonal, then the variations of the diagonalizing matrices V^'^^ do not contribute to the 
change in the eigenvalues and, according to first order perturbation theory 

5Ap(M) = {V^'hQ{M)V^'^^) 

= (y/°Vi5M[/2MoV/°)^)pp - (y/°ViMot/2Mo^5MMoV/°)^)pp (5.9) 
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If we introduce the matrices = 5MM^, v = v}^^ fiV}^^^ and Ki = v}^^UiV}^^^ and take 
into account that, due to ( ^.81) , 



we obtain the foUowing representation for ( |5.9| ): 

5Ap(M) = {K^uKlAiMo))^^ - (A(Mo)z.)^^ = - z.pp)Ap(Mo) (5.10) 



since A{Mq) is a diagonal matrix. 

The condition 5Ap{M) = means that {Kii'K\)pp = Upp for any p and any antiher- 
mitian u = —u^. This is possible only if Ki is diagonal, that is {Ki)pr = Cpdpr- Thus 
v}'^^ UiV}'^^^ = Ki where both Ui and Ki are diagonal matrices. This means that v}'^^ is 
a permutation matrix and Cp are the same as the eigenvalues {Ui)p (maybe in a different 
order). In turn, this implies the diagonality of UlVi'^^^ A{Mq)Vi'^^ and, because of ( |5.8|) , 
also the diagonality of M0U2MQ. Since U2 is diagonal, Mq also is a permutation matrix, 
and so is ¥2'^^ = Vi^'^Mq. In short, Vi^\ V^^^"^ and Mq are some permutation matrices. 

Going back to O we conclude that {vl^^UiVl^'^^){V^^'^U2V^^'^^)U3 = I. In other 
words, when properly ordered the eigenvalues of Ui, U2 and Us must give unity upon 
multiplication. 

This conclusion can be easily reformulated in terms of eigenvalue distributions ai{9), 
(72{0) and as{0). To do this, we introduce the "numbering functions" Sj(6'), j = 1,2,3, 
such that d'Ej{9)/d9 = CFj{6). These functions map the eigenvalue with number 
i G {1,...,A^} to the fraction G [0,1]. If we now construct the inverse mappings 
Vj{u) so that Ej(V^(u)) = u (no summation over j) then the multiplication condition for 
eigenvaluesH A^^^A^^'' Ap'* = 1 is equivalent to 



Vi{u) + V2{u)+Vs{u)=Q. (5.11) 



This is the selection condition imposed by the zero-area vertex for large A^. 

This discussion shows that the functions Vj [u) are somewhat analogous to momenta in 
a conventional field theory. Moreover, using the formalism of subsection 2.3 it is possible 
to show that the velocities Vj{9) in the collective field representation satisfy a similar 
constraint, 

vi{Vi{u)) =v2{V2{u))=vs{Vs{u)). (5.12) 



Here we implicitly presume that the eigenvalues have already been proprely ordered, so that 
in fact U1U2U3 = I is satisfied. 
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This constraint together with 



has a suggestive form. Let us recaU that the purpose of the zero-area vertices is to glue three 
cyhnders together. The functions Vj provide mappings between the boundaries of these 
cyhnders which are similar to the mappings between different coordinate patches forming 
the atlas of a manifold. This may mean that the string theory of two-dimensional QCD is 
most naturally formulated in terms of functions, functionally inverse to the collective field 
variables v and E. 

Another application of ( ^.11|) concerns the problem of the master field. The collective 



field theory gives us a way of evaluating the Wilson loops only without self-intersections. 
However, by joining a number of nonintersecting loops we can obtain an arbitrarily complex 
Wilson loop. Such joining is described by a vertex function we are discussing. However, 
to evaluate the composite Wilson loop one would need to know not merely the selection 
condition ( |5.11| ), but also the vertex function itself when this condition is satisfied^ This 
requires more subtle methods of calculation and is a separate subject that will be discussed 
elsewhere. 
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Appendix A. The change of variables in the collective field description of QCD^ 

Our goal is to find the large limit of equation ( pj.lll ). To this end we set Zn - 
exp N'^Fn and transform (|2.11|) into an equation for Fjy with the result 



where 



^dF^ 



2 



1 a . _ 1 ^ 



'k j^k 
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This would be a counterpart of the vertex coefficient in field theory. 



(A.l) 



Uj, = _^lnl)(^(i)) = _y cotan^ (A.2) 
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In the large limit all partial derivatives can be replaced by the derivatives with 
respect to the eigenvalue densities, 

OFn d SF 



N 



d0(}) d9 6ai{9) 



and all sums — by the integrals 



1 ^ 
N ^ 



27r 



fc=l 



ax{e) de. 



At the same time 



de 



(1)2 

k 



and can be neglected. As to the last term in (|A.1D, 



1 



N 



1 



d 



N 



N2 A^lQOH) 

k=l 



ATS 



1 



AT 



d 



N2 E ^^(1) 



(1)^ 



E 



AT 



J-'' 4 sin 



2 f^fc 



1 + l\-u 



k=l 



When N ^ oo 



and 



Uk ^ U{9) = 1 jcot^n(^-^^^ ai{9') d9' 



1 ^ 



27r 



A;=l 



t/2(6') ai{e) de. 



As to the term 



-y 



j^k 4 sm^ " , " 



it might appear that it is 0{1/N), since it contains only two summations (each of which 
naively contributes a factor of A^) divided by A^"^. However, due to the singularity at 
9k = Oj this term is in fact not negligible. Its 0{1) contribution comes from the regions 
where \k - j\ ^ N so that ol^^ - of^ ~ (A; - j) / {Nai{e\^^)) and 



-y 



ATS 



4 sin 

N 



2 "fe 



L V ' = J- f ff '^'"^C^") ^ 

^3 ,^(1) ^(l)xn /V3 



ATS 



j= 



27r 



fc = l n 
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where we used the identity 



Combining the pieces we obtain 

_ 27r 



dF r 
2^ = J M0)d9 



d 5F 



+ 2U{e) 



(A.3) 



If we now introduce a new functional S[ai{6), a2{0)\A^ according to (|2.12|) then 



d 5F 



+ U{9) 



d SS 



de6ai{e) ' " desaiie) 

so that (Q) entails (|2.13|) . 

Let us emphasize an important aspect of this derivation. In making the transition 
from ( A.l ) to ( |A.3p we relied on the fact that the large limit F = limAr^oo Fn exists. 
The same would not be true if we had attempted a shortcut, trying to replace V{6^^^)Zn 
by exp(A^^G'Ar) in (|2.11|) . In fact, the functional 'D{9^^^)Zn is not positive definite and 



therefore Gn does not have a well defined large limit. The functional 5" introduced 
2T2D is related to Zn by e^^{N^S) = \V{6^^^)\\V{e^^^)\ZN: the absolute value signs 



m 



being essential. In fact, it is the difference between T){6^^'^) and that gives rise to 



the interaction term ^ J a\{0) d6 in (|2.13|) 



Appendix B. The large limit of the Itzykson— Zuber integral. 



The large limit of the Itzykson-Zuber integral ( |2.32| ) can be studied using esentially 
the same technique that we used in appendix A to obtain the collective field theory of QCD. 
First, one represents In in the form 

. . X det||e^'"=''^'|| n\^n „2,y^iv j,2l det||e"^('"'-~^^')^|| 

In(A, B) = = e 2 Lz^fc=i "fe+Z^fc=i "^i !! 1. (B.l) 

Then it is easy to check directly that the quantity 

1 detlle"^^"''"^^)^ II 
A(afc)A(6j) 

satisfies the partial differential equation 

^^aj. 1 Aa^ 



dt A(a) ^ daf 

1=1 
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Once Jn is known, /^r can be retrieved as 



In = JN{t = 1) exp 



N 



N 



N 



k=l 



k=l 



Equation ( p.3| ) which is analogous to (|2.11| ) can be treated by methods described in ap- 
pendix A. The only difference is that the quantity Uk is given not by ( |A.2|) but by 



Uk 



1 d 
N'da'k 



In A 



with the large N limit 



U{a) 



a(a') da' 



a — a' 



(B.4) 



where a(a) is the density of the eigenvalues {afc}. Performing the transformations described 
in appendix A one deduces that the functional 



S\t\a,(3\= lim I — -r In JAr(t|A, S) l + - / a{a) a{a) \n\a ~ a'\dada 
N^oo j 2 J 



+ \ I l3{h)P{h') \n\h-h'\dhdh' 



(B.5) 



satisfies the following differential equation 



dS_ 
'di 



-\-oo 



-— = - a{a) 



d 5S 
da da{a) 



— a (a) 
3 ^ ^ 



da. 



(B.6) 



Like (|2.13|) this is the Hamilton-Jacobi equation for the dynamical system with the 
Hamiltonian 



H[p{a),U{a 



1 



+ 00 



p{a) 



da 



(B.7) 



However, as opposed to (|2.13| ), in this dynamical system p(a) is the distribution on an 
infinite real line rather than on a circle. In fact, it is the noncompactness of support of 
eigenvalue densities that distinguishes the Itzykson-Zuber integral from the large N QCD. 

Equation (|R6| ) is very hard, if not impossible, to solve in general. Fortunately, how- 
ever, it is easy to find its particular solution which describes the large N limit of the 
Itzykson-Zuber integral |T^]. This solution is given by the action of dynamical system 
( p.7| ) along its classical trajectory connecting the densities p{a) = a{a) and p{b) = (3{b) 
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within time t. It is not difficult to show that such quantity (which is known in classi- 
cal mechanics as the principal Hamilton function) does indeed satisfy ( B.6 ). Moreover, 
the variational derivatives at the end of the trajectory equal the corresponding canonical 
momenta: 

, 5S 

n(a,^ = o), — — = -n(6,t = 1). 



5a{a) ' ' " 6p{b) 
To find the above mentioned trajectory we must solve the equations of motion which follow 

from the Hamiltonian (|B.7|) . As in (|2.17|) , (|2.2g| ) these can be conveniently transformed 

>n(a 

da 



into a single equation for the quantity f{t, a) = -^^^ + inp{a): 



Then 

d dS _ dU{a,t = 0) _ _ ^ 

da 5a{a) da 



Now we are able to derive the constraints (p.36| ). We notice that the general solution 
of (|B.8|) can be written down in the parametric form 

x = RiO + FiOt 

B.9 

f{x,t) = F{0 

where -R(0 -^(0 some functions of the formal parameter ^. These functions should 
be determined from the initial conditions. The conditions to be imposed in our case are 
Im/(t = l,a) = 7ra(a) and Im/(t = 0, 6) = 7r/3(6). If we introduce the two analytic 
functions G+{x) and G-{x) according to 

G+(x) = x + f(t = 0,x) 

(B.IO) 

G- (x) = X — f{t = 1, x) 
then ImG'_|_(a;) = 7Ta{x) and ImG'_(a;) = —tt(3{x). In addition, from ( [B.9| ) we deduce 
G+{x) =x + [Fo R-^]{x) = [{F + R)o R-^] (x) 

as well as 

G-{x) = x- [F o {F + R)-^]{x) = [Ro{F + R)-^]{x) 
where o denotes functional composition. Therefore, 

G+{G-{x)) =G-{G+{x)) =x 

completing the proof of ( |2.36| ). 
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Appendix C. Wilson loops on a plane and the Hopf equation. 



In this appendix we will review the derivation of Wilson loops on a plane which uses 
the method of loop equations ||] |T^ . Using this example we will be able to see the close 
connection existing between the Hopf equation ( p.29|) and the Kazakov-Kostov equations 
for the two-dimensional QCD. In addition we will derive the explicit expressions ( 2.54 ) for 
Wilson loops on a plane. 



j / i j / i 




Fig. 6: The reconnection of a Wison loop at the self-intersection point employed 
in the two-dimensional loop equations. 



Generally, the loop equations relate the value of a self-intersecting Wilson loop W 
(fig. 7) to the values of Wilson loops with smaller number of self-intersections. In particular, 
if we "reconnect" the self-intersection on the left in fig. 6 to obtain the contour shown on 
the right then the Wilson loop W for the original contour satisfies the Kazakov-Kostov 
equation [§] 

RW = WW" (C.l) 

where the operator R is given byH 

f>=_^,A ^ ^ (C2) 

dSk dS, dSi dSj ' ^ ' ' 

In this formula 5'q, is the area of the window marked by the subscript a G {i^j, k, I}. 



We remember that in the two-dimensional QCD the Wilson loops depend only on the areas 
bounded by the contour but not on the contour shape. This follows from the invariance of QCD2 
with respect to the group of area preserving diffeomorphisms. 
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Fig. 7: A Wilson loop with two self-intersections (this corresponds to the winding 
number equal to 3). The areas of domains bounded by the loop are denoted as 
So, Si, S2, counting from the outside. 



Below we will be interested in the loops which wind several times around the circle of 
area A. In fig. 7 this would correspond to a situation when Sq, Si, . . . , Sn-i ^0, Sn ^ A. 
In this case the Wilson loop averages Wn{A) depend only on one variable A and, due to 
( |C.1| ), satisfy the equations 

(C.3) 



d 


d ^ 




dSo 






d 


d 






-iky- 


dSi 


dSo 




d 


d 


d 




dSk- 


-1 dSk+i 



( o o o \ 

( 2 = Wn-k-lWk 



These equations can be simplified if we introduce a new set of variables An = Sn, A^-i = 
An + Sn-i,...,Ak = Ak+i + Sk, . . .. Then d/dAo = d/dSo and d/dAk = d/dSk-d/dSk+i 
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so that from ( |C3| ) we obtain 



d 



d 



dAk dAk+i 



= Wr,-k-lWk, 



and, consequently, 



d d 



k-l 



dAo dAk 



) Vr„(Ao, Ai,...,A^) = J2 Wr,-k-i{A)Wk{A). (C.4) 



1=0 



Now, if S'o, 5"!, . . . , Sn-i then Aq = Ai — . . . — An = A and 



dA 



E 

k=0 



d 
dA, 



■Wn{Ao,Ai,...,A, 



(C.5) 



allAk=A 



Let us also recall that the dependence of on the external area Aq (or, the same. 
So) is always exponential, Wn ~ e~^°/^ so that 

oAq 2 



This, together with ( |C.4| ) and ( |C.5| ) allows us to write down a recursion relation for Wn- 

dWn n+1 



n-l 



dA 2 
Its solution has the general form 



(C.6) 



1=0 



Wn{A)=Pr,{A)e-'^'' 

where Pn{A) is a polynomial in A of degree n. It is easy to show by induction that these 
polynomials are given by ( 2.55| ) yielding the formula ( |2.54|) ||19|| . 

Our goal here is different. We would like to see how this formalism corresponds with 
the collective field description of QCD and, particularly, with the Hopf equation (|2.29| ). 

In fact, the recursion relation ( |(J.(j| ) resembles the Fourier transform of the Hopf equa- 
tion. Indeed, if we set 

oo 
n=0 



then ( |C6D entails 
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However, '0 is not the same as the coUective field function / because by construction is 
real while / has an imaginary part, Im/(t, 6) = na^{t, 9). But if we put C = —iO — ^ then 



Retl:(A,-ie - -) = ^Wr,{A) cos(n + 1)^ = 7ra*(t,^) - - 



n=0 



where we have used the relation between and 



27r 







(see (iH)). 

Thus the function 



has the correct imaginary part, Im/(A, 9) = Tra^{A, 9). Moreover, since 

the Hopf equation for t/j transforms into the Hopf equation for /: 

(^] +f(^] =0. 



(C.7) 



(C. 



Unlike '0 the function /, which has both real and imagnary parts, is precisely the function 
used in the collective field theory. Thus, the formula ( |C7| ) establishes the agreement 
between the Hopf equation and the loop equations of two-dimensional QCD. 
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